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We offer a new theoretical interpretation for the effect of enhanced electron density at '^Be nucleus 
encapsulated in fuUerene Ceo- Our ab initio Hartree-Fock calculations show that electron density 
at the ^Be nucleus in ^BeQCeo increase due to attractive effective potential well generated by the 
fuUerene. The 2s state in the isolated Be atom turns into 3s state in the joint potential. This 
new state has higher energy, and slightly larger amplitude at the Be nucleus than the previous 2s 
state. Moreover the 3s wave function has additional node appeared at the distance r ~ 5as from 
the center. The node imitates repulsion between the Be electron and the fuUerene wall, because 
the electron has zero probability to occupy this region. Such imitation of the repulsion by means of 
the node in attractive potential has direct physical analogy in the theory of a-a and A'^-A'^ nuclear 
interactions. 

PACS numbers: 23.40.-s, 21.10.Tg, 36.40.Cg, 61.48.+C 



I. INTRODUCTION 

In the ^Be /3-decay via electron capture (EC) process, the nucleus absorbs an electron from the atomic or molecular 
shell and is transformed to '^Li in the reaction p + e^ n + iy^- The decay rate is proportional to the electron density 
at the nucleus and therefore depends on the chemical environment of the radioactive isotope. ^Be has been used in 
investigations of the electron capture decay rate in various chemical states since the first studies by Segre [l| and 
Daudel [H . By now, experimentalists have published the results of more than sixty measurements relating to K- and 
L-shell electron capture by ^Be in different chemical forms and media. 

In 2005-2007, two teams studied the EC decay of ^Be inside the fuUerene Ceo Si'M- 1^ found, that the half-life 
of ^Be in metallic beryllium measured at room temperature exceeds the one in ^BeQCeo at room temperature by 
0.83% d, and by 1.5% Q if the latter is measured at 5°K. This difference between the ^Be EC /3-decay constants 
is the largest among available experiments. A density functional theory (DFT) based numerical calculations of the 
electron density at the Be nucleus have been presented in Q along with the experimental data. It was found that, 
in accordance with the experiment, the electron density at Be encapsulated in the center of Ceo is larger than the 
one at the nuclei in metallic beryllium. Qualitatively, in the metal, the electrons are shared, decreasing their local 
density at the nuclei, while the beryllium atom in Ceo remains intact. A more careful analysis, however, has shown 
that the electron density at ^Be@Ceo is larger by 0.17% even in comparison with that in an isolated beryllium atom. 
The authors explain this result as a "compression" of the Be's 2s orbital inside Ceo- The reason of the compression 
may be the "repulsive interaction" between Be and the Ceo cage according to [7j . 

II. EXAIVIINATION WITHIN THE FRAMEWORK OF HARTREE-FOCK BASED IVIETHODS 

Detailed understanding of the structure of atoms encapsulated in fullerens is important, in particular, for developing 
the concept of the fuUerene as an isolating cage which "does not affect" the trapped single atom and "protects" it from 
the outer environment. In the present work, we study the Be-Ceo interaction and its effect on the electron density at 
the Be nucleus within the framework of Hartree-Fock (HP) based methods. Although our value of the relative decay 
rate difference between metallic Be and BeQCgo is in qualitative agreement with the experimental one as well as with 
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the results of the previous theoretical studies, we suggest a new theoretical interpretation of the physical nature of 
the enhanced electron density at Be in Cgo- 

We started with the structural optimization of Be position inside Cgo- The fullerene's geometry was taken from 
experiment ^ (the length of the long and the short bonds are 1.448 A and 1.404 A, respectively) and fixed during the 
optimization, as the endohedral doping has, as expected, a small effect [7\. The total energy of BeQCeo complex at 
every trial configuration was calculated by the restricted (singlet spin state) Hartree-Fock method, with the 6-31G** 
molecular basis set [l3| in a cartesian form. Besides, the electronic correlations were taken into account within 
the second order perturbation method (MP2). For the calculations, we used our original program which employs 
the resolution of the identity (RI) method for the electron-electron interaction integrals and allows to perform the 
Hartree-Fock based calculations for large systems with moderate computational resources (see [ill, for details) . 
For both the Hatrtee-Fock and the MP2 variants, the optimization results in the position of the Be atom at the center 
of the fuUerene are in full agreement with the previous DFT based studies @, 0] • In order to evaluate the interaction 
energy defined as AE = E'soQCeo ~ (-E'bg + £-060 ) j have performed additional calculations of the total energies of 
Be atom, Ebc, and the fullerene, E'ceo- Besides, in a separate calculation, the basis set superposition error was taken 
into account by the counterpoise (CP) method [l^. The results are summarized in Table HI 



Table I: The Be-Ceo interaction energy in eV (Be is at the center of the fullerene). 





HF 


HF+MP2 


CP corrected 


0.91 


-0.41 


uncorrected 


1.00 


-0.63 



We therefore conclude that Be atom's equilibrium position at the center of the fullerene belongs to the attractive 
region of the Van-der-Waals interaction, i.e. the Be-Cgo interaction is attractive and the BeQCgo complex is stable 
in the ground state with respect to decay to Be and Cgo- This result is in contrast with the one of Lu et.al. Q 
who concluded a "slightly repulsive" Be-Ceo interaction from a DFT calculations, and obtained the value of -1-1.05 
eV for the interaction energy. We apply their speculations to our Hartree-Fock results, which also give the repulsive 
interaction as the pure HF method may not account of the dispersion energy. 
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Figure 1: Energy levels of molecular orbitals in Ceo, Be atom and BeQCeo- 



The consideration of Lu et.al. is based on the modifications in the energy levels of the HOMO/LUMO (highest 
occupied / lowest unoccupied molecular orbitals) of the Be atom and Ceo upon formation of the BeQCeo complex. 
Our results, presented in Figure [TJ differ qualitatively from those of FtJ, calculated with the B3LYP density functional 
with 6-31G** molecular basis set. In [21, Be's atomic 2s level lies in the middle of Ceo's LUMO-HOMO gap. Moreover, 
Ceo's LUMO-HOMO gap increases in Q upon endohedral doping with Be, because the HOMO is slightly lowered by 
0.03 eV, while the LUMO is elevated by 0.04 eV. These results were interpreted by Lu et al. as the evidence of a slight 
repulsion between the Be atom and the fullerene cage in the BeOCeo complex. Our calculations show that the Be's 
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2s orbital in the BeQCeo complex lies 0.14 eV lower than the Cgo's HOMO (see in Figure [T]). (The Be's 2s orbital in 
isolated Be atom lies 0.69 eV lower than the Cgo's HOMO.) Furthermore, Ceo's LUMO-HOMO gap does not change 
upon the endohedral doping, and both orbitals are just slightly lowered by 0.01 eV. Thus, there is little evidence 
of the "repulsive interaction", according to the definition of Lu at.al., between Be and Cgo within the Hartree-Fock 
calculations, even if the calculated interaction energy is positive. 

The Mulliken population analysis of our HF calculations with the 6-31G** molecular basis set gives the charge of 
+0.03 at Be compared with -0.14 in the DFT calculation of Lu et.al with the same basis set. This makes the concept 
of a slight hybridization of between Be and the fuUerene orbitals unacceptable within the Hatree-Fock method. 
However, strictly speaking, the Mulliken charges are not suitable basis set independent parameters for analysing the 
non-chemical interactions. 

Then, we consider the local electron density at the beryllium nucleus. The electron density for a polyatomic 
system is given by the standard relation, p(r) = I]„K„|(/3„(r)|^, where ipn{r) are normalized molecular orbitals, Hn ~ 
corresponding occupation numbers. In our case, the orbitals are the Hartree-Fock orbitals calculated with Danning's 
cc-pVTZ molecular basis set (Q) used in a Cartesian form. We employed two methods for the evaluation of the 
electron density at the nucleus which is formally located at the coordinate origin, r = 0. At first, we added narrow 
Gaussian s-functions (with exponents a between 10'^ and 10^ in inverse squared Bohr radius, a^) to the standard 
basis set. Then, we extrapolated the local electron density near the nucleus to the r — Q point by using the known (the 
consequence of the Kato theorem [l3|) non-relativistic asymptotic of the electron density for many-electron system at 
the Coulomb center with charge Z, {dp{r) / dr)r=o — — 2Zp(0). The coincidence of the both values of electron density 
at the Be nucleus with up to 0.02 a.u. (the atomic unit of density is equal to one electron in cubic Bohr radius, a^'^), 
i.e. about 0.05%, reflects the level of numerical error for our method. The calculated electron densities at single Be 
atom, metallic Be and BeQCeo are summarized in Table [III (Details of the procedures for calculations of electron 
density at metallic Be will be presented elsewhere.) 



Table II: Calculated electron density (in a.u.) at the Be nucleus. 

Orbitals 
1-st 2-nd Others Total 
BeOCeo 34.22 1.24 0.02 35.48 
Be atom 34.25 1.13 - 35.38 
Be metal 34.11 0.32 0.33 34.78 



Our results show that 
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i.e. a 2% decrease of the electron density at the nucleus from BeQCgo to metallic Be what is in qualitative agreement 
with the experimentally determined change of the decay rate at ^°K. Though, the absolute value is somewhat larger 
then the measured 1.5% as well as the value of 1.7% obtained by the DFT calculations [1, On the other hand 
our value is in excellent agreement with experimental data of Kraushaar et al. [l6| of direct measurement of the ^Be 

half-life in the metal source tIJ^ = 53.61 ±0.17 days and the half-life of ^Be in ^BeQCgo tI^^®^<'° = 52.47±0.04 
days from the work 

rp''Bc metal _ ji'^BciOCeo 

^ 7- ^ 100% ~ 2.1%. 

metal 

^1/2 

In the present paper we do not discuss these minor quantitative discrepancies neither between available experimental 
data, nor between various theoretical results. Experimental data require the further specification. As to numerical 
results it is necessary to take into account that both methods of calculation contain approximations. The Hartree- 
Fock method correctly takes account of the exchange, but lacks for the electronic correlations. The model density 
functionals, on the other hand, make approximations for both the exchange and the correlations. To make reasoning 
about the accuracy of the methods, especially in the case of the non-chemically, weakly bounded molecules, one would 
refer to a more robust theories like the coupled clusters or the configuration interaction methods. Instead, we are 
focussing on a qualitative phenomenon — the electron density difference between BeQCeo and isolated Be atom. The 
corresponding relative decrease of the electron density in our calculations is 0.28%, in qualitative accord with the 
value of 0.17% in [5, 6]. The reason of the enhanced electron density in BeQCeo might be a slight hybridization of 
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Be's and Ceo's orbitals, as suggested in [7|, so that the berylhum grabs the electron density from the fullerene. In that 
case, some of fullerene's orbitals would contribute to the electron density at Be. However, only two orbitals (those 
originated from atomic Is and 2s) make apparent contribution to the electron density at the Be nucleus. Thus, the 
hybridization concept has no support from the Hartree-Fock results. 

The structure of the Be is therefore changed due to the potential effect of the fullerene. In particular, fullerene's 
electrostatic field might prevent beryllium's electrons to spread out the cage acting as a strong repulsive potential 
wall. In that case, a more compact Is orbital would not be affected. However, the 2s one would rapidly vanish after 
certain distance from the center and would be compressed (due to the normalization) in the internal region resulting 
in the enhanced density at the nucleus. This concept, which also was expressed in Q, is in full accord with the data 
from Table Hn as well as with the increased energy of the 2s orbital of Be in Cgo- Nevertheless, a more detailed analysis 
of the spherically averaged electron density curves for the Is and 2s orbitals of atomic Be and Be in Cgo presented on 
Figure O rules out the concept of repulsive potential wall. Indeed, Be's 2s orbital inside the fullerene, though tends 
to zero at r ~ bas, then rapidly increasing and becomes even larger than the one of isolated Be atom. Below, we 
suggest a different solution to this intriguing issue. 
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Figure 2: Color online. Electron densities of Is and 2s states for isolated Be atom, and the first and the second orbitals of the 
Be atom in the BeQCeo complex: 1 (5) — 2s (Is) Be, Hartree-Fock; 2 (6) — 2s (Is) BeQCeo, Hartree-Fock; 3 (7) — 2s (Is) 
Be@C6o, model potential in Figure[3l 4 (8) — 2s (Is) BeQCso, model potential with model potential well in Figure|3] 

The electron density of the 2s orbital vanishes at r ~ Sas since the corresponding wave function crosses zero and 
changes its sign. Zero probability for the electrons to occupy the vicinity of r ~ 5as imitates the repulsive core for the 
2s electrons of Be in Ceo- However, the origin of the additional node in the wave function is the attractive potential 
well at 5aB < r < Sas . Spherically averaged electrostatic potential extracted from the Hartree-Fock electron density 
of Cgo (dashed-line curve in Figure [3|) is attractive. To illustrate the phenomenon, we designed a model spherical 
potential for the 2s orbital of Be (see in Figure[3]). It consists of the screened Coulomb potential as well as the spherical 
attractive potential well centered at r ~ 6.7aB- The screening constant was chosen to reproduce qualitatively the 
2s orbital for the isolated atom. It turns out, that the addition of the attractive potential into the model potential 
results in the appearance of the second node for 2s orbital, in the increase of its energy ( in full agreement with the 
result of HF calculation, see in Figure [T]) and also in the increase of the electron density at r = 0. 
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Figure 3: Color online. Potentials for the Be 2s orbital in the Be@C6o complex: 1 — Coulomb potential; 2 — model potential 
(screened Coulomb potential); 3 — potential well extracted from Hartree-Fock electron density; 4 — model potential "2" with 
model potential well. 

III. ANALYTICALLY SOLVABLE MODEL 

The reason of this phenomenon becomes clear from the following analytically solvable model. Let us compare the 
2s and 3s electron wave functions and the energy levels in the Coulomb potential (Figure H] (a) ) and in the new 
potential shown in Figure|4] (b). (This is the same Coulomb potential combined with the spherical potential layer.) 

We consider here only s wave functions ipns{r), because <pnp(0) = for all p states in non relativistic limit, and 
these p states do not give a contribution to the decay of ^Be. The electron radial wave functions for the combined 
potential in Figure 2] (b) are 

{aiexp(-Kr)iFi(l - Z/K;2;2Kr), < r < 
a2(sin(fcr) + &2 cos(fcr))/r, Ri < r < R2, 

as exp(-Kr)/r, R2 < r. 

Here, iFi is the confluent hypergeometric function, the wave numbers are k = y^2m\E\, k = \j2m(E — y), where 
E is the energy of the state (i5, V < 0). The coefficients ai_3, &2 and the energy E are determined from the conditions 
of continuity and differentiability of the wave function at the boundaries r — Ri and r = R2, where Ri,2 are internal 
and external radiuses of the potential layer correspondingly. 

For definiteness, we take charge Z = 2 for the Coulomb potential and depth V — —50 eV for the spherical layer 
located between i?i = 5.5 as and R2 = 7.5 as- The energy and the electron density at the nucleus of the 2s state 
in this potential are E2S — —13.61 eV and P2s(0) = 4 a^'^, correspondingly. If we add the spherical potential layer to 
the Coulomb potential, as it is shown in Figure S] (b) , and will gradually increase its depth up to the value of 50 eV, 
we will see the following. The energy of the 2s state decreases; its wave function gradually moves from the region of 
the Coulomb potential to the region of the spherical layer and becomes similar to the wave function of the Is state 
in the isolated potential layer shown in Figure |4] (c). The electron density of the 2s state at the nucleus decreases 
according to Figure [SI At the same time the 3s state gradually takes up the space region occupied previously by the 
2s state. Depending on the width of the spherical potential layer and on its depth V, the energy of the 3s state in 
the new combined potential may lie both above and below of the 2s energy level in the Coulomb potential, and may 
have either a higher or a smaller electron density at the nucleus (see in Figure [S]). 

Qualitatively, the mechanism of the variation of electron density of the 3s state at the nucleus, shown in Figure [SJ 
follows. The electron density in the origin, p„s(0), depends on the energy of the state Ens- For example, Pns(O) oc 
|£^„s|'^^^ in the Coulomb potential, Pns(O) oc \Ens\ in the infinite potential well and so on. The absolute value of the 
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Figure 4: Color online. Energy levels, electron wave functions, and electron densities of 2s and 3s states in: (a) Coulomb 
potential; (b) Coulomb potential combined with a spherical potential layer, (c) Energy level and electron wave function of Is 
state in a single spherical potential layer. 



energy, l-Ess], of the 3s state increases, when the depth \V\ of the spherical potential layer grows. At the same time, 
gradual redistribution of the 3s electronic density between the area of Coulomb potential and the spherical potential 
layer starts. For a relatively small \V\, the space occupied by the 3s wave function in the area of the potential spherical 
layer increases comparatively slowly, because the area of localization of the wave function practically does not vary, 
and the spherical potential layer does not have its own binding state. If the layer has a large depth, the situation 
changes. There is a binding state in the deep isolated spherical potential layer now. The energy of the 3s state in 
the joint potential verges towards the energy of such binding state, moving simultaneously from the position of the 
binding state in the pure Coulomb field. A part of the wave function occupies the forbidden for classical movement 
area (between the Coulomb potential and the spherical potential layer) , and the considerable enhancement of the 3s 
electron density occurs in the area of the spherical potential layer. The reduction of the wave function in the area of 
the Coulomb potential is not compensated anymore by growth of the |i?3s| energy, and the electronic density at the 
nucleus decreases fast, returning at first to the initial value (see in Figure [5]), and tends to zero at the further increase 
of the spherical potential layer depth. 

By applying the above considerations to the BeQCeo system one concludes the following. The Cgo fuUerene modifies 
the Coulomb potential of the Be atom in such way, that the 3s state in the joint potential occupies approximately the 
same position as the 2s state in the isolated Be atom. That is, the 3s state has approximately the same energy and 
practically the same electron density inside the Be atom as the 2s state in Coulomb potential. Moreover, corresponding 
wave function (p3s{r) has the second node, which imitates the repulsion of the electrons of the Be atom from the Ceo 
cage at r ~ Sa^. 

As regarding the Be 2s state, the energy of this state becomes considerably smaller in the new potential. Further- 
more, its wave function moves to the potential well formed by the fuUerene. 
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Figure 5: Color online. Electron densities of 2s and 3s states at the point r = as a function of spherical potential layer depth 
V. 



In that way, the small increase of electron density at the Be nucleus has casual character. One could obtain another 
result with different parameters of the "fuUerene potential well" (for example, by doping of Cgo with certain atoms). 
Thus, we obtain infrequent possibility to control the ^Be decay by means of a non-chemical interaction. In the 
considered particular case, this is the electrostatic interaction between the Be atom and the fuUerene. 

To confirm, that the proposed mechanism does not depend on the model, we also have considered two different 
analytically solvable models - a particle inside two "independent" and two connected spherical potential wells. The 
results obtained are in a good qualitative agreement with those described above. 



IV. COMPARISON WITH REPULSIVE CORE IN a-a AND N-N INTERACTIONS 



It is interesting to note in the end, that the problem considered here is not physically new and has a vague similarity 
with the well known problem of the repulsive core in nuclear physics. The concept of repulsive core in a-a and N-N 
interactions had been accepted as correct right until the seventies of the last century. This repulsive core arose at 
small distances as the consequence of the Pauli exclusion principle. It was established later that this concept was 
simplified. In some cases, the repulsive core must be understood in terms of the nodal wave function for relative 
motion in an attractive potential (see in Refs. and references therein). 

In other words, the zero probability for particle to occupy certain region can be achieved both by the infinite potential 
wall (repulsive core) and by the node of the wave function in an attractive potential. In this sense the phenomenon 
considered in the present paper resembles the above mentioned effects in the a-a and the N-N interactions. 



V. CONCLUSION 



In summary, according to our Hartree-Fock calculations with the electronic correlations accounted at the MP2 level, 
the lowest energy singlet configuration of Be atom encapsulated in Cgo is the one at the center of the fuUerene, - in 
accordance with previous DFT studied. The Be atom resides in the attractive region of the Van-der-Waals interaction 
with the interaction energy of about -0.6 eV, — in contrast with -1-1.05 eV from the DFT. Thus, the ^BeQCgo complex 
is stable with respect to decay to ^Be and Cgo- The HF electron density at the beryllium nucleus in BeQCgo exceeds 
the one in metallic Be by 2% in qualitative agreement with the relative difference of the corresponding ^Be EC decay 
rates measured by Ohtsuki et.al. as well as with the DFT calculations. The electron density at the Be nucleus in 
^Be@C6o also exceeds the one in isolated Be atom. The origin of this increasing is neither the modification of Be's 2s 
orbital in Ceo because of the hybridization with the fuUerene orbitals nor the repulsive potential wall at the region of 
the fullerene's atoms. Rather the replacement of the Be 2s state by 3s orbital in the new potential, which is the joint 
Coulomb potential of Be atom and the attractive effective potential well generated by the fuUerene. The 3s state has 
additional node at distance r ~ Sas from the center. This node imitates the repulsion between electrons of the Be 
atom and the Cgo cage, which has direct physical analogy in the theory of a-a and N-N nuclear interactions. 
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